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ABSTRACT 

Ultra high energy cosmic rays (UHECRs), with energies above ~ 6 x lO^^eV, 
seem to show a weak correlation with the distribution of matter relatively near 
to us in the universe. It has earlier been proposed that UHECRs could be ac- 
celerated in either the nucleus or the outer lobes of the nearby radio galaxy Cen 
A. We show that UHECR production at a spatially intermediate location about 
15 kpc northeast from the nucleus, where the jet emerging from the nucleus is 
observed to strike a large star-forming shell of g clS, IS cL plausible alternative. A 
relativistic jet is capable of accelerating lower-energy heavy seed CRs to UHECRs 
on timescales comparable to the time it takes the jet to pierce the large gaseous 
cloud. In this model many cosmic rays arising from a starburst, with a composi- 
tion enhanced in heavy elements near the knee region around PeV, are boosted to 
ultra-high energies by the relativistic shock of a newly oriented jet. This model 
matches the overall spectrum shown by the Auger data and also makes a pre- 
diction for the chemical composition as a function of particle energy. We thus 
predict an observable anisotropy in the composition at high energy in the sense 
that lighter nuclei should preferentially be seen toward the general direction of 
Cen A. Taking into consideration the magnetic field models for the Galactic disk 
and a Galactic magnetic wind, this scenario may resolve the discrepancy between 
HiRes and Auger results concerning the chemical composition of UHECRs. 
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Introduction 



To understand the origin of cosmic rays (CRs), it is important to distinguish between 
the lower energy CRs which can be contained within the magnetic field of our Galaxy and 
thus have energies of up to about 3 x 10^^ eV for heavy nuclei, and those that are even 
more energetic. The bulk of the CRs below that energy can be explained by supernova 
explosions, while the extr emely energetic ones probably originate from either some c lass 
of Active Galactic Nuclei ( iGinzburg &: Syrovatskiilll964t iBiermann &: Strittmatterl 119871 ) or 
some extreme type of stellar activity such as gamma-ray bursts ( 1Waxmanlll995l ) . Indeed, the 
spectrum of CRs shows a kink near 3 x 10^® eV, matching the expectation that their origin 
changes around this energy threshold. 

Stellar explosions can account for the flux, spectrum, particle energy and chemical 
composition of the less ener getic CRs, cons i dering that all very massive stars exp lode into 



jxpl( 

their pre-existing winds (e.g. lPrantzoslll984l : IStanev et al.lll993t iMeyer et al.lll997l ). Further 



quantitative confirmation of this picture has now emerged from detailed observations of 



cosmic ray electrons and positrons, as well as the WMAP haze ( IBiermann et al.l l2009a 



2010l ). The supernova ori gin of Galactic cosmic ra ys may lead us to an understanding of the 
seed particle population ( iBiermann et al.ll2009bl ) on which active galactic nuclei energizing 
radio galaxies can operate their acceleration processes. 

The origin of ultra high energy cosmic rays (UHECRs) is still an unresolved issue, but 
a few clues have begun to emerge. Although their arrival directions are nearly isotropic, a 
general correlation with the distribution of matter has been noted by the Auger observatory 



(IStanev et al. 



1995 



tory ( lAbbasi et al 



2008 



Auger-CoU.I l2008al Jbl). although it is disputed by the HiRes observa- 



2010a|). In particular, there may be excess events with arrival 



direc tions close to the nearby radio galaxy Centaurus A ( lAuger-Coll.ll2008at [Abraham et al. 



20091 ). There are contradic ting claims from expe riments as to whether the U HECR events 



are he avy nuclei (Auger - [Abraham et al.l |2010| ) or purely protons (HiRes 
2010bl ). Both possibilities need to be explored. 



Abbasi et al. 



In a picture where UHECR eii ergies are attained by a single kick up from a seed popu- 
lation ( IGallant fc Achterberg|ll999l) through the ac tion of a relativistic jet, these events can 
indeed involve heavy nuclei (IBiermann et al.l 2009bl) . In such a scheme the seed particles are 



the CRs near the spectral knee (IStanev et al.lll993[ ) and the relativistic shock is very likely to 
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arise from a jet carving out a new channel after being launched from a primary central black 
hole that has been reoriented followin g the merger of the nuclear black holes of two merging 
galaxies (iGergely fc Biermannll2009f). In this scenario all the UHECR particles are a mix of 
heavy nuclei, and the spectrum in IStanev et al.l (Il993[ ) actually gives an adequate fit to the 
Auger data ( Biermann et al. 2009bl). See Fig. 1. T he sky distribution is easily isotropized 
by the intergalactic magnetic fields (jPas et al.ll2008l ): f or the case of heavy nuc lei one is even 
confronted with the possibility of excessive scattering (iBiermann et al.l 2009bl ). This picture 
also allows the incorporation of the Poynting flux limit ( iLovelacd 119761 ): the particles to be 
accelerated must remain confined within the jet diameter. This condition translates into 
a lower limit for the jet power, allowing most UHECR particles to originate from the jet 
interacting with lower energy CRs produced in the starburst in the central region of Cen A. 

Therefore we explore a scenario bas ed on the observed head-on enco u nter of the Cen A 



jet with magnetized interstel lar clouds (IGopal-Krishna &: Saripallil Il984l : iKraft et al.l 12009 



Gopal-Krishna fc Wiitall2010l ) from which UHECR acceleration ensues. A distinctly appeal- 
ing aspect of this proposal is that the postulated jet-cloud interaction is actually observed 
within the northern lobe of Cen A, whereby the jet is se e n to be disturbed, bent westward and 
possibly disrupted temporarily ( iMorganti et al.l Il999l : lOosterloo &: Morgantil l2005l ). Since 
any supersonic flow reacts to a disturbance with shock formation, this in turn could cause 
particle a cceleration. Note that the Fermi/LAT error circle for the peak of the gamma-ray 
emission (lAbdo et al.ll2010l ) encompasses the jet-cloud interaction region, at the base of the 
northern middle lobe of Cen A, about 15 kpc from the nucleus. 



2. Acceleration in Cen A from a jet interacting with gaseous shells 



The key point is that the interaction of the northern jet with a gaseous shell in the 



northern middle lobe has clearly been seen (lOosterloo fc Morgantil l2005l : iKraft et al.ll2009l ) 



and massiv e star formation is revealed at the location of the interaction by the GALEX 
UV image (IKraft et al.ll2009l ). Although other mechanisms can bend and disrupt radio jets, 
onl y a jet-shell interaction can explain the variety of data (radio, HI, UV, X-rays) for Cen 
A ( iKraft et al.ll2009l ). It has also been argued that the oft-debated peculiar morphology 
of the no rthern middle radio lobe can be readily understood in terms of the same jet-shell 
colhsion ( IGopal-Krishna fc Wiitall2010l ). 



An important aspect of the basic acceleration physics to be stressed is that when parti- 
cles are accelerated in a shock pro pagating par a llel to the magnetic field lines, the maximurn 
particle energy Emax is given by ( iHillasI Il984l : iGinzburg fc Syrovatskiil Il964l : IStanevI |2004| ) 
Emax = G Z PshRs B, where e is the elementary electric charge, Z is the numerical charge 
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of the particle, l3sh is the shock speed in units of the speed of hght, the available length scale 
is Rb, and the strength of the magnetic field is B. However, when the shock propagation is 
highly oblique, the corresponding limit (jJokipiilllQSTt iMeli fc Biermannll2006l ) becomes 



cZRbB, 



(1) 



which is independent of the shock velocity. Invok ing relativistic shocks obvio usly adds an 
additional factor of 'jsh, the shock's Lore ntz factor (IGallant &: Achterberelll999l). Loss es will 
curtail this maximum attainable energy (jHillaslll984t iBiermann &: Strittmatterl 119871 ) . 



We now focus on the particle acceleration due to the observed interaction of the jet 
with shells of fairly dense gas. Gen A has long been known to have a nu mber of stel 



lar sh ells, located in the vicinity of both the Northern and Southern lobes ( iMalin et al. 
19831 ). Some of these shells hay e later been found t o contain large amounts of dense 



atomic (jSchiminovich. et al.lll994( ) and even molecular fjCharmandaris et al.l |2000| ) gas (~ 
7.5 X 10^ M0). These shells are gen erally though t to have originated from the merger of a 
massive elliptical with a disk galaxy fjQuinrull984| ). very probably the same merger that gave 
rise to the peculiar overall appearance of this large elliptical galaxy marked by a striking 
dust lane. Radio maps reveal that the northern jet has encountered such shells at dis- 
tances of 3.5 and 15 kpc from the core, and fiared up eac h time to the same side, thereby 



forming the northern-inner an d the northern- middle lobes (IGopal-Krishna fc Saripallilll984 



Gopal-Krishna &: Wiita]l201Cll ). Simulations of such c ollisions indicate the formation of strong 



shocks where the jets impinge upon gas clouds (e.g. iGhoi et al.l 120071 ) 



We must ask whether the maximum observed particle energies, of order 10^^ eV, are 
actually attainable in such interactions. Accelerating particles to such copious energies 
requires that the Larmor motion of a particle must fit within the gaseous cloud, both before 
and after the shock that forms inside the cloud by interaction with the impingir ig relativisti c 
jet. This leads to the condition Emax ~ e Z Rd, also called the Hillas limit ( lHillaslll984l ). 
which is a general requirement to produce UHEGR via shocks. 

Adopting the very reasonable parameter values of 3 kpc for Rd , the approximate ob- 
served size of the HI shell found in th e Northern Middle Lobe of Gen A (jOosterloo fc Morganti 
2OO5I : iGopal-Krishna fc Wiital[201ol ). and 3 X 10 ^ Gauss for the magnetic field, it follows 
that the ene rgy must remain below Z x lO^^eV. Since particles are observed up to about 
3 X 10^" eV flBird et al.lll994l ). this implies that heavy nuclei, such as Fe, are much preferred 
for this mechanism to suffice; however, if a stronger magnetic field were present, this would 
ease the requirement on the abundances and allow for GRs to be accelerated to even higher 
energies. The magnetic field in the shell is not well constrained, but the required value is 
modest. 
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The Hillas limi t condition f lHillasI Il984( ) mentioned above can be expressed another 
way (jLovelacd Il976l ). Taking the energy needed for particle acceleration to derive from a 
jet, we can connect the time-averaged energy flow along the jet with the condition that the 
accelerated particles are contained within the jet diameter, 



>10"^erg s- 



int 



E„ 



Z X 102ieV 



(2) 



where fint is an intermittency factor describing the temporal fluctuations of the energy 
outflow. Equality in this critical expression would imply that the energy flow in the jet is 
an entirely electromagnetic Poynting flux, an unrealistic extreme scenario. For Cen A we 
require both an intermittency factor < 1, and presumably also heavy nuclei, e.g., Z 26. We 
find /int ~ 0.75 in order to match the kinetic jet p ower, which has been argued to be Lj^t 
10^'^ e rg s~^ through se veral different approaches (jWhysong fc Antonuccil l2003l : lAbdo et al. 



2010 



Kraft et al.l 



(lAharonian et al 



2009 



, and references therein). The recent HESS observations of Cen A 
20091 ) detected an ultra-high energy (> 250 GeV) photon luminosity of 



2.6 X 10^^ erg s ^, but the entire photon luminosity in gamma-rays (> 100 keV) is 



only 
~ 2 X 10^2 grg s 



and thus also consistent with L 



•jet 



10^3 



erg s 



We next examine whether the inferred luminosity of UHECRs is indeed attainable. 
Assuming the observed spectrum of the jet corresponds to a CR particle spectrum of about 
E~'^''^, this leads to the requirement that the observed power in UHECR particles must be 
multiplied by a factor of about 200 in order to integrate over the power-law spectrum. The 
data then require a luminosity of about 10^^ erg s~\ still below the inferred jet power of 10'^'^ 
erg s~^ for Cen A f Whysong fc Antonucci 2003 : Abdo et al.lboiol ). Thus, we could allow 
for a duty cycle of 0.1, and still have adequate jet power. So the jet's interactions with a 
dense cloud are capable of powering the observed UHECRs. Another way of asking the same 
question is, can a jet actually catch a sufficient number of particles from the knee region 
with energies near PeV and accelerate them to the ankle region near EeV to ZeV? Assuming 
that the energy density of CRs in the starburst region is about 100 times what we have in 
our galaxy, the particle density near and above 10^^ eV is about 10~^^ per cc. If through the 
non-steadiness of the jet these CRs are caught at the same rate by a kpc scale jet having 
an opening angle of, say, 5°, the cross-section of ~ 10^^'^ cm^ implies a rate of 10"^^ particles 
accelerated per second. Pushing them to UHECR energies gives an energy turnover of order 
]^g42 gj,g g-i j^g^ ^YiQ energies above 10^*'^ eV, again quite sufficient. 

Third, we need to check whether enough time is available for the particles t o be accel- 
erate d. A jet encounter with such a large cloud would last for at least 10^ yrs (jChoi et al. 



20071 ). A shock in either the external or the internal medium would take some small multi- 
ple of the Larmor time scale at the maximum energy of a few times 10^'^ yrs, to complete 
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the acceleration process. The two relevant time scales, for transit and acceleration, seem 
consistent within the scope of our broad estimates. 

Lastly, we need to check whether the time scales are long enough so that the time 
window for possible detection of the UHECR source is not too brief. The time scales for 
particle acceleration and the jet-cloud encounter are somewhere between lO'^ and 10^ yr. 
The times for the jet to transit a shell and then to move on to the next shell appear to be in 
a ratio of about 1 to 10. Therefore, a duty cycle, fint, of about 0.1, which is easily allowed 
for by the above calculation, is actually necessary to maintain a quasi- continuous output of 
accelerated particles. 



3. Consequences of the jet/cloud acceleration scenario 



Having shown that the basic model is viable, we now consider some of its consequences. 

First we no te that the jet may st i ll be mildly precessing after the episode of the merger 
of black holes f lGergely fc BiermannI |2009| ) , i n the afterm ath of the merger of the ellip- 
tical and spiral galaxies comprising Cen A ([Israel Il998l ). Also, the gaseous shell may 
have its own motion, also due to the preceding merger of the two galaxies. This would 
na turally explain the observed multiple bendings and flarings of the northern jet in Ceri 
A JCopal-Krishna &: Saripallilll984l : ICopal-Krishna et al.lliooa ICopal-Krishna fc Wiitalboioh . 
Both effects would expose continuously fresh material to the action of the jet, but are not 
an essential requirement for our model. 

Second, the transport and scattering of the particles along the way might smooth out 
any variability even if Cen A were the only significant source of UHECRs in our part of 
the uni verse. Such variability might explain the i nconsistencies between Auger and HiRes 
results (lAbraham et al.ll2010l : lAbbasi et al.ll2010al ). The magnetic field at the site of origin 
is locally enhaiiced by the Lorentz factor of the shock, possibly between 10 and 50 (e.g., 
Biermann et al.l l2009bl ) . That could imply the shortest possible variability time r^ar — 
100 Ti,ar,2 y^s, taking a high Lorentz factor of 50. The scattering near the Earth needed to 
attain near isotropy in arrival directions requires a relatively strong magnetic field within the 
distance equal to cr^ar- So the containment of Fe particles of up to 3 x 10^° eV would imply 
an energy content near Earth of Es^var ~ 6 x 10^^ r^^2 ^rg. Interestingly, this total energy 
approaches the energy of a hypernova (10^^ erg). However, there is currently no evidence for 
such a region surrounding the Sun. 



Finally, we have to follow through with the deduction from the Poynting flux limit 
(jLovelacelll976l ). that the highest energy events can only be heavy nuclei if they come from 
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Cen A. This limit requires that all particles caught by a shock in the jet have E/Z less than 
or equal to that of Fe at lO^"'^ eV, the highest energy event yet seen; let us assume initially, 
that this one event at 10^°'^ eV is a factor of 3 below the real limit imposed by the acceleration 
site, the shock in the jet interaction region. It follows that He above 10^^'^, and CNO above 
]^q2o.i gY are ruled out, but near 10^^ '' eV both are possible. We use the prescriptions of 
AUard et al.l ( l2008l ) to define a photo-disintegration distance Kdis for any nucleus and energy. 



Averaging over some wiggles in the curves that cover both FIR and microwave backgrounds 
from the very early universe, we find that over the relevant energy and nucleus charge range 
an adequate approximation is h^dis = 10^'^Mpc {Z/Zpe)/{E/10^^''^ eVY'^ , which we use here 
to guide us. There are two extreme scattering limits. In one limit, the isotropization of 
the events from Cen A is done iii thejntervening intergalactic medium (IGM). Cosmological 
MHD simulations by iRyu et al.l ( 120081 ) imply a Kolmogorov approximation; then Atrav = 
10^-^ Mpc [{Z / Zpg) / {E /lO^^'"^ eV)Y^^. However, this already leads to extreme losses of the 
heavy nuclei between Cen A and us. So we consider the other limit, in which the UHECRs 
travel es sentially straight fr om Cen A to us, and are isotropized in the magnetic wind of our 
Galaxy ( lEverett et al.ll2008l ). Modeled values of the wind's magnetic field strength (~ 8 /iG) 
and radial scale (~ 3 kpc) allow Fe, as well as all elements down to about Oxygen, to be 
scattered into isotropy; however, there is less effect on lower Z elements. No other approach 
gave a reasonable fit to the data. The losses due to the path traversed during the scattering 
are small. A fit with this approach is shown in Fig. 1. One c ould use other magnetic wind 
model numbers, closer to a Parker-type wind (jParkerl Il958l ). but the essential results do 



not change. Now we must ask, how ca n this be compatible with IGM models (IRyu et al. 



2008; 


Das et al. 


2008; 


Cho & Rvu 


2009) 



IGM, s cattering can be reduced if much of the overall magnetic energy is pushed into thin 



sheets ( iBiermann et al.l l2009d ) and such substructure plausibly arises from radio galaxies 
and galactic winds. A second question is whether the magnetic field could also produce a 
systematic shift on the sky for UHE CRs, in add i tion t o scattering and isotropizing them. 



Indeed, any Galact ic magnetic f ield (iBeck et al.l 



of a Galactic wind ( IStanevlll997l : Everett et al.l 120081 ) would also produce a systematic shift 



19961 ). in the disk or in the foot region 



relative to the central position of Cen A on the sky. Since Cen A is not far from the sensitivity 
edge of the Auger array in the sky, it is quite possibl e that there is a s h ift fo r all events, 
especially at slightly lower energies. The models of IZirakashvili et al.l ( 119961 ) show that 
angular momentum conservation and transport quickly generate a magnetic field component 
parallel to the galactic disk, which would shift particle orbits in a direction perpendicular to 
the disk, and possibly away from the center of symmetry. 

A testable prediction of this scenario then is that a solid angle on the sky containing 
half the UHECR events towards Cen A should show the signature of lighter nuclei, hence 
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larger fluctuations, compared to the events seen from the remaining part of the sky. Since the 
main scattering also has a systematic component, the center of this anisotropy may be shifted 
with respect to Cen A, so that the part of the sky with the largest fluctuations in the shower 
properties may be offset by up to a few tens of degrees from Cen A for Z > 1. This effect 
might be strong enough so that in some parts of the sky lighter elements might predominate 



over heavies and thus reconc ile results from the Auger and HiRes experiments (lAbbasi et al. 



201 Obi : [Abraham et al.ll2010l ) . This could soon be checked with the growing data on UHECRs. 
If such a test were positive, it would unequivocally and simultaneously show that Cen A is 
the best source candidate, that scattering depends on the energy/charge ratio, and that the 
most energetic events are heavy nuclei. 
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Fig. 1. — Testing the shift of the spe ctrum in paper (|Stanev et aLlll993l ) with Auger data, 
using the propagation calculations of ( lAUard et all 120081 ) with a 3.8 M pc distance to Cen A 
and assuming all isotropizing is in the magnetic wind of our Galaxy (lEverett et al.ll2008l ). 
Note that the break here is due to the MHD structure of massive star winds, pushed to EeV 
energies by a highly relativistic shock. A color version appears in the on-line journal. 



